Temperature and pressure-induced spin-state transitions in LaCoOs 
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We report the continuous variation of the spin moment of cobalt in LaCoOs across its temperature 
and pressure-induced spin transitions evidenced with K/3 emission spectra. The first thermal tran- 
sition is best described by a transition to an orbitally nondegenerate intermediate spin (5=1) state. 
In parallel, continuous redistribution of the 3d electrons is also indicated by partial fluorescence 
yield X-ray absorption spectra. At high pressure, our study confirms that the material becomes low 
spin between 40 and 70 kbar at room temperature. 

PACS numbers: 71.30.+h,71.70.Ch,75.30.Wx,78.70.En 



I. INTRODUCTION 

The perovskite-type cobalt oxide LaCoC>3 has at- 
tracted much attention since the 1950s due to its ther- 
mally induced spin-state and insulator-metal transitions. 
LaCoOa has a non-magnetic insulating ground state, 
with all the six 3d electrons of Co 3+ occupying the ti g 
level because the crystal field splitting exceeds the in- 
traatomic exchange energy. Consequently, Hund's rule 
does not apply, and the cobalt is in low-spin (LS) state 
with S—Q for the t\ g ^ g configuration. At higher tem- 
peratures two broad transitions have been observed in 
the magnetization measurements, occurring at around 
90 K and 500 K. The classic interpretations invoked 
the appearance of the t% g e 2 g (S = 2) high-spin (HS) state 
of Co 3+ at the first transition, and various ordered su- 
perstructures of the HS and LS cobalt atoms at higher 
temperatures!^ Inconsistencies started to arise, how- 
ever, when neutron diffraction studies failed to find any 
superstructures, 3 and more conflicting discrepancies were 
brought up by experimental studies of the electronic 
structure. Namely, it was found that neither Co 2p and 
O Is X-ray absorption spectroscopy (XAS) nor photoc- 
mission spectra show pronounced variations below 300 K; 
they reflect a gradual transition only above 500 K from 
the LS state to a mixed spin state with a main contri- 
bution from the HS stated The failed observation of the 
90-K transition and the conflicting description of the tem- 
perature behavior stimulated significant further efforts 
to determine the magnetic and spin states of LaCoOs. 
Shortly after these experiments, the 90-K transition was 
evidenced by a spin-polarized neutron-scattering study. 5 
Not much later, a completely different description was 
proposed that included the formation of a t\ g e g (•5=1) 
intermediate-spin (IS) state at the first transition** Theo- 
retical works completed this LS /IS scenario with propos- 



ing an orbital ordering for the IS state, which melts at 
the second transition£& This has become the prevail- 
ing scenario, and most of the recent experimental obser- 
vations are interpreted as supporting the IS stateA±2iiI 
Although including this state in the explanations seems 
to explain the observed phenomena, some authors still 
invoke the HS state at higher temperatures 
Another sign of the incompleteness of a simple LS /IS sce- 
nario is that the description of the high temperature sus- 
ceptibility data usually needs a treatment different from 
the low temperature ones, e.g. invoking change in the 
orbital degeneracy at the second transition^ Thus, the 
interpretations of the two transitions are still conflicting; 
the precise spin state as a function of temperature, and 
the variations of the electronic structure associated with 
the transitions are not yet fully understood, and require 
further studies. 

Pressure can also cause spin-state transitions in tran- 
sition metal compounds^ as the crystal field splitting 
rapidly increases with the shrinking of the bond lengths 
(10 Dq oc r~ 5 ). Since the transitions of LaCoOs are 
accompanied by anomalous volume expansions, it is ap- 
parent that relatively low pressure can have considerable 
effect on the spin state. The effect of pressure to the 
90-K transition was studied up to 18 kbar by Asai et 
at, who found that the energy gap between the LS and 
the higher spin state increases with pressured Recently, 
Vogt et al. performed an X-ray powder diffraction study 
to examine the pressure- induced changes. The results 
were interpreted as a continuous transition from IS to LS 
state, which is completed at around 40 kbariiS Chemical 
pressure, introduced by partial substitution of La 3+ with 
the smaller Eu 3+ , leads to similar stabilization of the LS 
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We present here high-resolution X-ray spectroscopy 
study of the spin state and the electronic structure of 
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cobalt in LaCo03 as a function of temperature and pres- 
sure. The applied techniques allowed us to follow the 
occupancy and the local magnetic moment of the cobalt 
3d orbitals. The applicability of the X-ray absorption 
near-edge structure (XANES) spectroscopy to study the 
electronic structure is well known, as it gives access to 
the unoccupied density of states. In the case of the spin 
transition, it is expected that the pre-edge of the main 
Is — > 4p absorption line, which has a relevant contri- 
bution from quadrupolar Is — > 3d transitions, follows 
sensitively the variations in the populations of the ti g 
and e g levelsiSi These orbital populations are expected 
to undergo relevant changes, as the spin transition in 
a localized picture is best described as redistribution of 
the electrons between the tig and e g orbitals, as it was 
described above. The applied partial fluorescence yield 
(PFY) detection technique^ even if it is not a true mea- 
sure of the absorption cross section^ allows better en- 
ergy resolution than conventional XANES,— and can be 
used to monitor more sensitively the population of the 
different statesj 2 ^ 2 ^ X-ray emission spectroscopy (XES), 
in fact, the Kj3 emission line can be utilized to probe the 
local magnetic moment on the 3d orbitals, i.e., the lo- 
calized spin moment i^iiSSiSi The origin of the sensitivity 
to the spin state is the strong exchange interaction be- 
tween the 3p hole and the 3d orbitals in the final state of 
the emission process; the two different possible orienta- 
tions of the spin of the 3p hole with respect to the resul- 
tant local moment of the unpaired 3d electrons leads to a 
splitting. 28 The arising spectra are not simple due to the 
numerous multiplet terms involved; nevertheless, a low- 
energy satellite (K/3 1 ) to the main emission line (K(3\^) 
is observed, whose intensity is proportional to the spin 
localized on the 3d orbitals^ Besides this phenomeno- 
logical description, a full theoretical treatment based on 
crystal-field or charge-transfer multiplet theory can sat- 
isfactorily reproduce the spectral variations with the spin 
state^ii 



II. EXPERIMENT 

LaCoC>3 was synthesized as described in Ref. ^3- For 
reference samples, high purity LiCo02 and C0F3 were 
purchased from Sigma- Aldrich. The experiment was car- 
ried out at the ID16 beamline of the European Syn- 
chrotron Radiation Facility. The undulator radiation 
was monochromatized with a cryogenically cooled Si(lll) 
double-crystal monochromator, and focused to a spot size 
of 150x100 /im 2 with a Rh-coated toroidal mirror. Ar- 
ranged on a Rowland-circle spectrometer of 1 m diame- 
ter, a spherically bent Si(620) crystal and a Peltier-cooled 
Si diode were applied to analyze and detect the emitted 
radiation, respectively. The K(3 spectra were collected by 
scanning the Bragg angle of the analyzer around 70.69°, 
which corresponds to the maximum of the Co K(3 emis- 
sion line; the detector was also translated and rotated ac- 
cordingly. The overall resolution of the spectrometer was 



1.5 eV. Fluorescence yield XANES spectra were taken 
with a photon-counting Rontec detector. Partial fluo- 
rescence yield XANES spectra were taken on the max- 
imum of the _fC/3-emission line by scanning the incident 
energy through the absorption edge (and keeping the ana- 
lyzed energy constant)^ For the temperature dependent 
studies, a scattering angle of 90° was chosen. The low- 
temperature experiment was carried out in a standard 
closed-cycle helium cryostat. For the high-temperature 
studies, the sample was mounted on a Ta/Mo support 
heated by a ceramic heater device, in a vacuum cell. 

A Mao-Bell diamond anvil cell with a Re gasket was 
used for the high-pressure experiment. The sample was 
prepared by grinding a piece of a LaCoC>3 crystal, and 
a methanol-ethanol (4:1) pressure medium was added 
to the powder sample to ensure hydrostatic pressure. 
The pressure was determined by the ruby fluorescence 
method. Both the incident and the scattered radiation 
passed through the diamond anvils; the scattering an- 
gle was 30°. In order to minimize loss of the emitted 
radiation, a diamond of smaller height was used on the 
downstream part of the cell: the diamonds had heights 
of 2.5 mm and 1.5 mm, with culet sizes of 0.7 mm and 
0.6 mm, respectively. The energy of the exciting beam 
was set to 13 keV, a compromise between the absorption 
loss in diamond and scattering cross section. 

All spectra were normalized to the beam monitor; the 
XES spectra were further normalized to the spectral area, 
and the spectra were aligned to have the main lines at 
the same energy. 



III. RESULTS 

Figure ^ shows the XES spectra of our HS (C0F3) 
and a LS (LiCoG^) reference compounds. As it was ex- 
pected from previous observations jSLSii and theoretical 
calculations^ the largest spectral variations are seen on 
the intensity of the main Kf3\^ line and on its K/3' low- 
energy satellite (see the difference spectrum in Fig.^). 

Similar variations, though to a lesser extent, were ob- 
served in the line shapes of the spectra of LaCoC>3 (Fig. 
|2J), as a function of temperature. At the lowest tem- 
peratures the spectra are identical to that of LiCoC>2. 
As temperature increases, the increase of the satellite 
intensity and the decrease of the main line intensity is 
observed. These variations are continuous with temper- 
ature, which reflects a continuous change of the average 
spin of the cobalt. The 913-K spectrum exhibits a sudden 
change. While below this temperature the spectral vari- 
ations were reversible, the last spectrum remained the 
same when the sample was cooled down after collecting 
data at 913 K. This indicates that a chemical transforma- 
tion took place when reaching the highest temperature of 
our study; this will be discussed later in more detail. 

The PFY XANES, displayed on Fig. shows also con- 
tinuous variation with temperature. Below the Is — > Ap 
main absorption line at 7726 eV, several low intensity fea- 
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FIG. 1: K/3 spectra of LiCo0 2 and CoF 3 , a LS and a HS 
Co 3+ compound, and their difference. 
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FIG. 2: (color online) Evolution of the LaCoOs K/3 spec- 
tra with temperature, arrows indicate the main line shape 
variations with increasing temperature. At the bottom, dif- 
ferences from the lowest temperature spectrum are shown for 
17 K, 150 K, 300 K, 503 K, 803 K, and 913 K. 

tures appear. Following earlier interpretations) 3 ^* 3 ^ 3 , the 
broad feature at around 7716 eV has been assigned to a 
ligand-to-metal charge transfer (LMCT) shake-down pro- 
cess, which results in a final state of ls 1 3d 7 4p 1 L, instead 
of the ls 1 3<i 6 4p 1 of the main transition, (L denotes a hole 
on the ligand). At even lower energies, the pre-edge re- 
gion reflects transitions to the ti g and e g levels, partly 
through direct quadrupolar transitions, but also through 
to dipolar transitions to d-p hybridized orbitals. The lat- 
ter will affect the pre-edge intensity but not the spectral 
shape, which still shows characteristic features of the 3d 
crystal- field splitting even in the dipolar component^* At 
higher temperatures, the edge progressively shifts while 
the pre-edge region, as it will be discussed later in more 
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FIG. 3: (color online) PFY XANES of the LaCoOs on the 
max of K/3, with the pre-edge region in the inset. Arrows 
indicate the spectral variations at increasing temperature. 



detail, reveals a continuous redistribution of electrons on 
the tig and e g levels through intensity changes, such as 
broadening and the shift of their spectral weight to lower 
energies. At 913 K, similarly to the XES, a drastic change 
was observed in the spectrum: the absorption edge was 
shifted by several eVs to lower energies, and the pre-edge 
was also shifted and reduced. The large edge shift indi- 
cates a change in the valence state of cobalt, which was 
confirmed by chemical analysis: the sample was trans- 
formed to LaCoC>2.5, which reveals that a Co 3+ — ► Co 2+ 
transformation and simultaneous oxygen loss took place. 
This is not entirely surprising, as Radaelli and Cheong 
have found a third lattice anomaly (besides the two at 
the magnetic transitions) around 800 K, which was at- 
tributed to formation of oxygen vacancies. 14 Moreover, 
such reduction is known to happen to LaCoOs at even 
lower temperatures in the presence of a reducing agent. 35 
As only a few spectra were taken at low tempera- 
tures, rapid total fluorescence yield XANES scans were 
measured between 15 K and 300 K to verify that the 
spectral changes (hence the redistribution of the 3d elec- 
trons) are indeed continuous in this region. The spectra 
are shown in Fig. ^ the lowest and highest-temperature 
spectra are similar to those of Toulemonde (taken at sim- 
ilar temperatures)) 3 ^ in between a continuous change is 
seen as displayed in the inset of the figure. The pre-edge 
was interpreted as a result of two contributions: 32,33 a 
low-energy feature of small, but increasing intensity was 
attributed to transitions to the ti g orbitals, while the di- 
minishing larger feature at higher energy was attributed 
to e g orbitals. The relative change of their spectral 
weight is therefore a sign of the ti g — > e g electron trans- 
fer. Although this is a simplified description, the nature 
of the changes of the electronic structure indeed originate 
from the redistribution of the 3d electrons; therefore, for 
monitoring purposes this description should suffice. As 
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FIG. 4: TFY XANES at 15 K and 300K. The intensity vari- 
ation of the low-energy (i,2 g ) shoulder in the pre-edge of the 
Co K-edge is shown in the inset. Above this plot, a spectrum 
taken at T=300K illustrates (after background subtraction 
and normalization) the fitting with which the plotted values 
are obtained. 
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FIG. 5: Evolution of the LaCoOs K/3 spectra with pressure. 
The main figure shows the extremes and their difference, while 
the inset displays the K0 satellite region for all the measured 
pressure points, which include 0, 7, 22, 42, 70, and 100 kbar. 



IV. DISCUSSION 



A. Model calculations 



we see in the inset of Fig. 0J which presents simple two- 
Gaussian fits to the pre-edges, the variation of the fea- 
tures of the pre-edge are indeed continuous. To derive 
a quantitative description of the redistribution of the 3d 
electrons, a detailed resonant X-ray emission study would 
be desirable. 



As for the pressure dependence, a gradual variation of 
the K(3 line shape was found up to 70 kbar. As it is seen 
in Fig. |SJ the relative intensity of the satellite were re- 
duced while that of the main line were increased, which 
is typical for a transition toward the LS state. Above 
70 kbar no further changes were observed. In this case 
we could not monitor the redistribution of the 3d elec- 
trons directly with XANES due to the strong absorption 
of the diamonds at the Co K edge. However, the nature 
and the extent of the changes in the XE spectra between 
ambient pressure and 70 kbar are the same as those which 
are observed when the temperature is decreased at am- 
bient pressure from 300 K to 17 K; this similarity will be 
more obvious in section llV Bl where this comparison can 
be made in a more quantitative manner. At 70 kbar and 
above the spectra appear the same as the lowest temper- 
ature one of the temperature dependent series. This sug- 
gest that at room temperature a pressure-induced tran- 
sition to the LS state is essentially completed between 
40 kbar and 70 kbar, which is in agreement with the 
findings of Vogt et ali& 



In order to facilitate the understanding of the variation 
of the experimental data and to help to relate it to the 
spin state, the Co Kj3 emission and the pre-edge at the 
K-edge of Co in LaCoOs were modeled using a charge 
transfer multiplet approach 3 ^!. The XANES spectra 
were also calculated using an ab initio multiple scattering 
code. 

The charge transfer multiplet model used a 3c? 6 + 
3d 7 L + 3d s L 2 ground state basis set, where L denotes a 
hole in the ligand valence band. The Slater integrals were 
obtained by the Hartree-Fock method 38 and scaled by 0.8 
to take into account intra-atomic configuration interac- 
tion effects. For LS couplings atomic values were used. 
The crystal field perturbation in Oh symmetry was con- 
sidered using the approach developed by Butleo 3 ^, and 
charge-transfer effects were taken into account by using 
the code by Thole and Ogasawara™. The Is lifetime 
was set to 1.3 eV and the final state lifetime broadening 
FWHM was approximated as 0.2 x Ef, where Ef is the 
final state energy^i. 

Although the IS state cannot be obtained as a true 
ground state in the multiplet approach, it was assumed 
that using the triplet 3 Ti term of 3d 6 can provide a rea- 
sonable approximation for the spectra of this state. The 
model parameters for the different spin states of LaCoC>3, 
as well as for a high-spin and a low-spin reference com- 
pound of 3d 6 configuration, CoF 3 and LiCoC>3, are pre- 
sented in Table [I] The parameter values were obtained 
from the literature^iMii^ii^ 6 . and they were modified 
to some extent when necessary to achieve better agree- 
ment with the experimental data. The calculation with 
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TABLE I: parameters of the multiplet calculations. For de- 
scription of parameters see the text. 
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FIG. 6: Calculated Kf3 spectral shapes of LaCoOa for the 
different possible spin states; the mean of the HS and LS 
spectra is also shown for comparison with that to the IS. 
The differences from the LS spectrum are also plotted in the 
bottom of the figure. 



the LS ground state was compared to the experimen- 
tal spectrum of LaCoOs measured at 17 K and to that 
of LiCo02- The calculation on the intermediate spin 
ground state was tentatively compared with the experi- 
mental spectrum measured at 403 K, and C0F3 was used 
to represent also LaCoOa in the high-spin state. In the 
table A stands for the charge transfer energy, defined as 
the energy difference between the centers of gravity of 
d n and d n+1 L configurations and Udd is the Hubbard U 
given by E(d n+1 ) + Eid 11 ' 1 ) - 2E(d n ). The core hole 
potential Udc is set to the same value for both Is and 
3p holes. V e is the hybridization strength for the e g or- 
bitals in the ground state; the hybridization strength for 
t2 g states was set to half of the V e value. For states with 
a core hole the V e hybridization strength was reduced by 
0.4 eV. 

The simulated Kf3 spectra for the different spin states 
of LaCoC>3 are shown on Fig. The spectra are shifted 
so that the main lines match the nominal Co Kf3 en- 
ergy. In agreement with the experimental data and with 
previous observations^^ and calculations^ the inten- 
sity of the K(3' low-energy satellite of the main K/3i ^ 
line is increasing with the increasing cobalt spin moment. 
However, there are no apparent spectral variations that 
would provide an easy way to verify the presence of the IS 
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FIG. 7: (color online) Co K edge of LaCo03 calculated with 
FEFF8.2 using structural data from literature^- An arrow in- 
dicates the variations at increasing temperature. 



state. According to Tsutsumi's proposal, the position of 
the satellite with respect to the main peak should move 
with the spin state as AE = J(2S + 1) (where J is the 
3p — 3d exchange integral)^ thus for the IS a substan- 
tial upward shift (almost 6 eV) of the K f}' feature would 
be expected. However, as observed also at other spin 
transitions, 21 this is not the case: no new features show 
up, and the calculated IS spectrum barely differs from 
the mean of those of the HS and LS (c/. Fig. [JjJ. This 
is also in accord with the temperature dependent exper- 
imental data on LaCoOs provided that the first transi- 
tion is indeed an LS— >IS one, since no relevant energy 
shifts of the satellite were observed. In general, the vari- 
ations of the experimental spectra (when available) are 
well accounted for by the calculations. The calculated 
line shapes suffer from only a slight overestimated height 
of the a Kj3' satellite, which originates from the approx- 
imations of the lifetime broadening model and the over- 
estimation of transition matrix element strengths in the 
calculation. However, due to the discussed problems, the 
satisfactory agreement is not sufficient to extract quan- 
titative information about the spin state or discriminate 
between LS/HS or LS/IS scenarii in LaCoOs. We will 
show that careful line shape analysis can overcome this 
shortcoming (c/. section HV B|l . 

For the description of the Co XANES spectra, an ab 
initio multiple scattering approach was used to study 
the spectral variations with the temperature for a few 
selected values between 10 K and 700 K. The temper- 
ature dependent structural data was taken from the 
literature^. The program utilized was FEFF8.2^£i^ a 
real-space multiple scattering code. For the calculations 
the Hedin-Lunqvist self-energy and a cluster of 134 atoms 
was used. Dipolar and quadrupolar transitions were both 
considered to get a reasonable estimation for the pre- 



6 




7.708 7.71 7.712 7.714 7.716 7.718 
Energy (keV) 



FIG. 8: Results of the multiplet calculations on the pre-edge 
of the Co K-edge of LaCoOs for the different spin states 

edge. As seen from the simulated spectra shown on Fig. 
the calculations reasonably describe the edge shifts, 
which accompany the structural changes with tempera- 
ture. However, as expected for a one-electron approach, 
it fails to reproduce the LMCT peak, and does not show 
relevant variation on the pre-edge. Therefore, for a bet- 
ter description of the pre-edge region a different theoret- 
ical approach, the charge transfer multiplet model was 
invoked. The quadrupole contributions were calculated 
using the multiplet code with the same parameters as pre- 
sented in Table [I] The dipolar part was modeled using 
a phcnomenological density of p symmetry final states 
similar to the model used by Caliebe et al^ in calcu- 
lating the pre-edge of Fe K- absorption edge of Fe2C>3. 
The charge transfer multiplet calculation predicts a sin- 
gle peak in the pre-edge for the LS state, which originates 
from transitions into e g levels. When higher spin states 
are populated, the lower-lying t-x g levels also become ac- 
cessible at the absorption process due to their partial 
depopulation. The appearing features are expected to 
cause a 'widening' of the pre-edge, especially in the di- 
rection toward lower energies. This is in agreement with 
the experimental data presented in Figs. |3]andQJ The 
determination of the spin transition on the basis of this 
spectra is again challenged by the relatively small differ- 
ences of the IS and HS line shapes. 

B. Line shape analysis 

Despite the conspicuous effects of the localized spin 
on the Kj3 spectra, and the evident origin of the satel- 
lite, the task of determining the spin from the line shape 
is not obvious. As multiplet theory shows, the spectra 
are made up by contributions of numerous terms, and 
the arising lines can even have different widths due to 
the mentioned different term-dependent lifetimes. There- 



fore, fitting the spectra with theory would lead to many 
parameters whose number might be close to that of the 
data points. For a phenomenological description, the Kj3 
spectra can usually be fitted reasonably with three Voigt 
functions^ This might even serve as a useful approach 
for examining trends over a large set of compounds, for 
example; nevertheless, it has little use in quantifying the 
spin state for a transition which is not completely known. 
A more effective approach is to use the combination of 
proper reference spectra of known spin states of tran- 
sition metals with similar local structure to model the 
spectral shape, or eventually to fit it to the data. In such 
cases the spectra can be reasonably reproduced, and the 
extent of the transition can be procured, as it was demon- 
strated in Ref. lU However, this requires that the refer- 
ence spectra is taken under similar conditions, i.e., simi- 
lar resolution, source size, background, and it is seldom 
that these are all satisfied. Moreover, it might not always 
be easy to find appropriate reference materials, e.g., for 
the present study no IS Co 3+ compound was available, 
and using the CoF 3 as HS reference can also give less 
satisfactory results for a Co06 system due to the differ- 
ences in the bonding, let alone the uncertainties faced 
when interpolating the line shapes for the not entirely 
known IS state. So far the most successful approach is to 
follow the variation of the spin state with this technique 
was to plot the integrals of the absolute values of the dif- 
ference spectra {hereafter IAD), and then scale them to 
values deduced from some known reference materials. 51 
The IAD value for spectra SP^ with respect to a reference 
spectra SP re / is obtained as numerical approximation to 
Eq.HJ 



IAD ( = / | SP 4 - SP re/ | dE, (1) 

(Note that this integral would give without taking 
the absolute values due to the normalization of the spec- 
tra.) This approach deserves further consideration as it 
has many advantages. It is model independent, it can be 
performed without invoking results of theoretical calcu- 
lations, it relies only on proper reference spectra or IAD 
values obtained for a similar transition, and the raw data 
needs only a simple preparation before the analysis can 
be used, so it is a promising tool for the rapid determi- 
nation of the spin state from the XES spectra. On the 
other hand, it requires that the whole set of experimental 
data are acquired under similar conditions. 

The area differences were so far used for the determi- 
nation of the spin state in a pragmatic manner^i Nev- 
ertheless, it can be easily demonstrated that in case of a 
transition between two states, this approach, although it 
might look arbitrary at first, is a true probe of the line 
shape variation, which is parallel to the variation of the 
spimSi 

In the present case, however, three spin states might 
be involved in the transitions. In order to follow the spin 
state through the line shape variations, we have to exam- 
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FIG. 9: IAD values for the calculated and the measured 
reference spectra. 
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ine, whether and how the IAD correlates with the spin. 
Unfortunately, we have no means of expressing this in 
form of analytical functions, we can only rely on values 
obtained from the simulated and the experimental line 
shapes. In Fig. [§] the experimental IAD values for the 
A5=0, 1, 3/2, and 2, as well as the IADs derived from 
the calculated spectra of Fig. for the S=0, 1, and 2 
states are presented. The IAD obtained for the C0F3 HS 
reference compound (compared to the LiCo02 LS refer- 
ence) is 0.12 for the AS=2; LaCo02.5, which is known to 
have S=3/2, gives an IAD of 0.084. Finally, for a molec- 
ular spin-transition system of Co 2+ with a A5=l tran- 
sition an IAD of 0.049 was obtained^ Apparently both 
IAD sets vary linearly with the spin, thus on the basis 
of the existing data it seems to be at least a promising 
approximation to determine the spin state. The linear- 
ity possibly arises from the the similarity of the environ- 
ments of Co, as the ligands provide a somewhat distorted 
octahedral surroundings in all cases, and it is well known 
in multiplet theory that the local symmetry is more im- 
portant to determine the spectral shape than the nature 
of ligands or the valence stated The larger values of the 
calculated ones can be explained by the overestimation 
of the Kfi' intensity mentioned above; rescaling this data 
to the experimental one using the fraction obtained from 
the proportion of the slopes of the fitted lines, suggests 
that the calculation overestimates the line shape varia- 
tion by 17%. 

We followed this approach to analyze the data; the 
IAD values obtained with Eq. ^ are presented in Fig. 
IT01 In order to convert these values to spin, the linear 
dependence derived from Fig. were used to rescale the 
data to determine the average local spin of cobalt as a 
function of temperature or pressure. The obtained spin 
values are presented on the right scale of Fig. Note 
that the estimated spin of LaCo02.s ( a t T=913 K) is 
consistent with its known value. 



FIG. 10: Integrals of the absolute values of the area differ- 
ences of K[3 spectra. Top: temperature dependence; the open 
and filled symbols represent data from two different samples. 
Bottom: pressure dependence at room temperature. The dot- 
ted lines represent one half increments in the spin value. 

C. Testing the data with the different scenarios 

To determine which spin-transition scenario is in agree- 
ment with the data, the temperature behavior of mag- 
netic susceptibility or anomalous volume expansion are 
usually analyzed in terms of thermally excited localized 
two or three-state modelsii2*i^2i££ In terms of the three- 
state model, we can express the average spin on cobalt as 
the sum of the contribution of the (non-zero spin) excited 
states: 

2 

S(T) = J2 S l v,{2S l + l)e- A */ fcB7 >, (2) 

i=l 

where Si is the spin of the excited state i, is its 
energy with respect to the LS ground state, and i>i is its 
orbital degeneracy, while Z denotes the partition func- 
tion. Thus, in the LS/IS/HS scenario Si = 1,S 2 = 2, 
while Vi can have different values; for the two-state mod- 
els the unnecessary term is omitted. 

Equation |21 was fitted to the obtained spin data (pre- 
sented in the top of Fig.llO[) to determine the model and 
orbital degeneracies compatible with it. In Fig. ^2 we 
compare our experimental results on the spin of cobalt 
and fits to the various models. The simple two-state mod- 
els do not give an appropriate description, especially the 
fits with the LS/HS scenario, and those of the LS/IS sce- 
nario with higher IS orbital degeneracy (uis > 1) were 
so poor that they are not included in the discussion. As 
even the remaining possibile two-state model, which in- 
volves a transition to the IS excited state with i/ig = 1, 
gives a poor fit at high temperatures (see the dashed 
line), a model with a change in the IS orbital degeneracy 
(i^is = 1 — > 3) at high temperatures was also tested. Such 
"switching on" of the orbital degeneracy at high tern- 
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FIG. 11: (color online) Fit of the temperature evaluation 
of the cobalt spin state with two- and three-state models de- 
scribed in the text. The legend describes the different orbital 
degeneracies assumed in the particular models. 



TABLE II: Model orbital degeneracies and fit parameters for 
the spin on cobalt. The first line describes the LS/IS two- 
state model, the second is the LS/IS scenario with a change 
in the orbital degeneracy above 400K. The other lines belong 
to LS/IS/HS three-state models. 

vis ms Ais / meV A H s / meV Aig / K A H s / K 



1 




16.2 




188 




1- 


^3 - 


16.6^32 




193^371 




1 


1 


16.7 


130 


194 


1504 


1 


3 


16.6 


180 


193 


2085 


3 


1 


31.8 


187 


369 


2164 


3 


3 


31.8 


223 


369 


2592 


6 


1 


43.2 


253 


501 


2938 


6 


3 


43.2 


278 


501 


3223 



perature was suggested by Zobel et al. to attain a bet- 
ter agreement with susceptibility data^S However, this 
model, which is shown in Fig. ^2 as a dashdot line with 
a jump at 400 K, does not bring a satisfactory improve- 
ment to the fit of the high temperature data. Among 
the three-state-models, the ones with an orbitally non- 
degenerate (z^is = 1) IS first excited state follow most 
reasonably the temperature evolution of the spin state 
(solid lines). The necessity of invoking the third state is 
evident from a comparison of these lines with that of the 
previously described two-state models. The fits with the 
other three-state models that involve higher orbital de- 
generacy increase too rapidly, so they give poor descrip- 
tion for the entire temperature region (dotted and dashed 
lines). With these models, varying the orbital degeneracy 
of the HS state does not yield a relevant difference, and 
thus, for clarity, only those with (fis = 3), which give a 
slightly better fit, are included in the figure. In brief, 
in contrast to recent susceptibility measurements ^SiSi 
the data is best described with the three-state model. 
The parameters of the fits are shown in Table |Hj the de- 
scription of its lines follow the description of the lines in 
Fig. ^2 presented above. In summary, three-state mod- 
els (LS/IS/HS) with an orbitally nondegenerate (^=f) 
IS first excited state give the best description of the tem- 
perature evolution of the spin state. The energy gaps 
obtained for the v=\ IS cases, which correspond to tem- 
peratures of 188-194 K, are in good agreement with the 
values of 180 K and 185 K, obtained by Zobel et al. re- 
cently from magnetic susceptibility measurements and 
from anomalous thermal expansion, respectively^ and 
with the 180 K obtained with 59 Co NMRpS 

It should be noted, however, that we assumed the en- 
ergy gaps Ai to be independent of temperature. Albeit 
this conservative approach was applied in many recent 



works with success ^2d2£l modifications of the crystal- 
field splitting and hybridization strength with tempera- 
ture are likely to affect the energy of the excited states. 
Hence, based on a description proposed long ago^ Asai 
et al. applied corrections to the spin-state energies, which 
varied with the extent of the transformation^ while 
Baier et al. suggested a straightforward way to extract 
the temperature dependence from the experimental data 
at the cost of restraining themselves to the two-state 
model. In the latter, the extracted A(T) exhibits a 
moderate variation at low temperatures, while above the 
insulator-metal transition (T>500 K) it drops rapidly. 20 
Applied to our data with the appropriate reformulation of 
Eq. the results of this approach suggest different high- 
temperature behavior (namely, a more rapid decrease of 
A(T) above room temperature) than that obtained from 
the susceptibility datai^i Although we do not favor an- 
alyzing the data with this approach due to its strong 
restrictions, we also refrain here from using the other 
description^ since temperature-dependent gaps would 
require introducing more fit parameters, and this would 
be challenged by the limitations of the data set presented 
here. Nevertheless, with the "conservative" approach 
a reasonable description of the experimental data was 
achieved, thus we expect that the variation of A, should 
not be dramatic. Actually, many aspects of the temper- 
ature behavior of LaCoC>3 support this (especially below 
the insulator-metal transition), which include relatively 
small variations of the relevant bond lengths^ small ex- 
cess entropy^ and relatively small changes in the vi- 
brational spectra^ 5 . A more detailed study on LaCoC>3 
and its partially substituted compounds will follow this 
work, which should permit us to address the nature of 
the temperature dependence of the higher spin state(s) 
using more refined models. 
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V. CONCLUSIONS 

We have studied the temperature and pressure induced 
spin-state transitions of LaCoC>3 with high resolution x- 
ray spectroscopies. The XE spectra reflects continuous 
temperature variation of the spin localized on the cobalt. 
Similarly, the PFY XANES spectra shows the contin- 
uous redistribution of the 3d electrons between the t2 g 
and e g levels of cobalt. We have developed a model- 
independent way to quantitatively follow the spin tran- 
sition with X-ray emission spectroscopy through the in- 
tegrated absolute difference spectra. Our high-pressure 
study confirmed that the pressure-induced transition re- 
ported by Vogt et al. 19 is indeed a transition to the LS 
state. The temperature behavior of the cobalt spin was 
found to be most compatible with a transition to an or- 
bitally nondegenerate IS state with a gap of about 190 K. 



The higher transition requires invoking new parameters 
into the description; our data was best described with a 
partial transition to a HS state. Future detailed experi- 
ments utilizing resonant and non-resonant X-ray emission 
spectroscopy should shed more light on the variations of 
the electronic and spin state. 



Acknowledgments 

The authors thank Prof. L. H. Tjeng, Prof. M. Abd- 
Elmeguid and Dr. G. Monaco for stimulating discussions, 
Drs. T. Lorenz and C. Zobel for providing the LaCo03 
samples, Dr. S. Huotari for his careful reading of the 
manuscript, Dr. R. Verbeni for the preparation of the an- 
alyzer crystal, and C. Henriquet for technical assistance. 



1 J. B. Goodenough, J. Appl. Phys. 29, 513 (1958). 

2 P. M. Raccah and J. B. Goodenough, Phys. Rev. 155, 932 
(1967). 

3 G. Thornton, B. C. Tofield, and A. W. Hewat, J. Solid 
State Chem. 61, 301 (1986). 

4 M. Abbate, J. C. Fuggle, A. Fujimori, L. H. Tjeng, C. T. 
Chen, R. Potze, G. A. Sawatzky, H. Eisaki, and S. Uchida, 
Phys. Rev. B 47, 16124 (1993). 

5 K. Asai, O. Yokokura, N. Nishimori, H. Chou, J. M. Tran- 
quada, G. Shirane, S. Higuchi, Y. Okajima, and K. Kohn, 
Phys. Rev. B 50, 3025 (1994). 

6 R. H. Potze, G. A. Sawatzky, and M. Abbate, Phys. Rev. 
B 51, 11501 (1995). 

7 M. A. Korotin, S. Y. Ezhov, I. V. Solovyev, V. I. Anisimov, 
D. I. Khomskii, and G. A. Sawatzky, Phys. Rev. B 54, 5309 

(1996) . 

8 P. Ravindran, H. Fjellvag, A. Kjekshus, P. Blaha, 
K. Schwarz, and J. Luitz, J. Appl. Phys. 91, 291 (2002). 

9 S. Noguchi, S. Kawamata, K. Okuda, H. Nojiri, and 
M. Motokawa, Phys. Rev. B 66, 094404 (2002). 

10 C. Zobel, M. Kriener, D. Bruns, J. Baier, M. Griininger, 
T. Lorenz, P. Reutler, and A. Revcolevschi, Phys. Rev. B 
66, 020402(R) (2002). 

11 G. Maris, Y. Ren, V. Volotchaev, C. Zobel, T. Lorenz, and 
T. T. M. Palstra, Phys. Rev. B 67, 224423 (2003). 

12 K. Asai, A. Yoneda, O. Yokokura, J. M. Tranquada, 
G. Shirane, and K. Kohn, J. Phys. Soc. Jpn. 67, 290 
(1998). 

13 D. Louca, J. L. Sarrao, J. D. Thompson, H. Roder, and 

G. H. Kwei, Phys. Rev. B 60, 10378 (1999). 

14 P. G. Radaelli and S.-W. Cheong, Phys. Rev. B 66, 094408 
(2002). 

15 C. N. R. Rao, M. M. Seikh, and C. Narayana, Top. Curr. 
Chem. 234, 1 (2004). 

16 T. Kyomen, Y. Asaka, and M. Itoh, Phys. Rev. B 71, 
024418 (2005). 

17 M. P. Pasternak, R. D. Taylor, R. Jeanloz, X. Li, J. H. 
Nguyen, and C. A. McCammon, Phys. Rev. Lett. 79, 5046 

(1997) . 

18 K. Asai, O. Yokokura, M. Suzuki, T. Naka, T. Matsumoto, 

H. Takahashi, N. Mori, and K. Kohn, J. Phys. Soc. Jpn. 



66, 967 (1997). 

19 T. Vogt, J. A. Hriljac, N. C. Hyatt, and P. Woodward, 
Phys. Rev. B 67, 140401(R) (2003). 

20 J. Baier, S. Jodlauk, M. Kriener, A. Reichl, C. Zobel, 
H. Kierspel, A. Freimuth, and T. Lorenz, Phys. Rev. B 
71, 014443 (2005). 

21 G. Vanko, T. Neisius, F. Renz, S. Karpati, 
A. Shukla, A. Mirone, and F. de Groot, 

ESRF Highlights pp. 59-60 (2002), URL 

http : //www . esrf . f r/UsersAndScience/Publications/Highlights/2( 

22 K. Hamalainen, D. P. Siddons, J. B. Hastings, and L. E. 
Berman, Phys. Rev. Lett. 67, 2850 (1991). 

23 P. Carra, M. Fabrizio, and B. T. Thole, Phys. Rev. Lett. 
74, 3700 (1995). 

24 P. Glatzel and U. Bergmann, Coord. Chem. Rev. 249, 65 
(2005). 

25 C. Dallera, M. Grioni, A. Shukla, G. Vanko, J. L. Sarrao, 
J. P. Rueff, and D. L. Cox, Phys. Rev. Lett. 88, 196403 
(2002). 

26 J.-P. Rueff, C.-C. Kao, V. V. Struzhkin, J. Badro, J. Shu, 
R. J. Hemley, and H. K. Mao, Phys. Rev. Lett. 82, 3284 
(1999). 

27 J. Badro, J.-P. Rueff, G. Vanko, G. Monaco, G. Fiquet, 
and F. Guyot, Science 305, 383 (2004). 

28 K. Tsutsumi, H. Nakamori, and K. Ichikawa, Phys. Rev. 
B 13, 929 (1976). 

29 G. Peng, F. de Groot, K. Hamalainen, J. Moore, X. Wang, 
M. Grush, J. Hastings, D. Siddons, W. Amstrong, 
O. Mullins, et al., J. Am. Chem. Soc. 116, 2914 (1994). 

30 X. Wang, F. M. F. de Groot, and S. P. Cramer, Phys. Rev. 
B 56, 4553 (1997). 

31 M. G. Kim, Y. S. Im, E. J. Oh, K. H. Kim, and C. H. Yo, 
Phys. B 229, 338 (1997). 

32 O. Toulemonde, N. N'Guyen, F. Studer, and A. Traverse, 
J. Sol. Stat. Chem. 158, 208 (2001). 

33 O. Haas, R. P. W. J. Struis, and J. M. McBreen, J. Sol. 
Stat. Chem. 177, 1000 (2004). 

34 M.-A. Arrio, S. Rossano, C. Brouder, L. Galoisy, and 
G. Calas, Europhys. Lett. 51, 454 (2000). 

35 B. G. Tilset, H. Fjellvag, and A. Kjekshus, J. Sol. Stat. 
Chem. 119, 271 (1995). 



10 



36 F. M. F. de Groot, A. Fontaine, C. C. Kao, and M. Krisch, 

J. Phys. Condens. Matter 6, 6857 (1994). 51 

37 F. M. F. de Groot, Chem. Rev. 101, 1779 (2001). 

38 R. D. Cowan, The Theory of Atomic Structure and Spectra 
(University of California Press, Berkeley, 1981). 

39 P. H. Butler, Point Group Symmetry Applications: Meth- 
ods and Tables (Plenum Press, New York, 1981). 53 

40 See the following website, URL 
http : //www . anorg . chem.uu.nl/people/ staf f/FrankdeGroot/" 4 

41 M. Taguchi, T. Uozumi, and A. Kotani, J. Phys. Soc. Jpn. * 
66, 247 (1997). 

42 T. Saitoh, T. Mizokawa, A. Fujimori, M. Abbate, 56 
Y. Takeda, and M. Takano, Phys. Rev. B 55, 4257 (1997). 

43 M. Abbate, R. Potze, G. A. Sawatzky, and A. Fujimori, 
Phys. Rev. B 49, 7210 (1994). 

44 A. Chainani, M. Mathew, and D. D. Sarma, Phys. Rev. B 
46, 9976 (1992). 

45 T. Matsushita, M. Mizumaki, N. Ikeda, M. Nakazawa, 
A. Agui, Y. Saitoh, T. Nakatani, A. Yoshigoe, and S. Naka- 
mura, Surface Review and Letters 9, 1327 (2002). 

46 T. Mizokawa and A. Fujimori, Phys. Rev. B 53, R4201 
(1996). 

47 A. L. Ankudinov, S. D. Conradson, J. Mustre de Leon, and 
J. J. Rehr, Phys. Rev. B 57, 7518 (1998). 

48 A. L. Ankudinov, C. E. Bouldin, J. J. Rehr, J. Sims, and 
H. Hung, Phys. Rev. B 65, 104107 (2002). 

49 W. A. Caliebe, C.-C. Kao, J. B. Hastings, M. Taguchi, 
A. Kotani, T. Uozumi, and F. M. F. de Groot, Phys. Rev. 
B 58, 13452 (1998). 

50 S. D. Gamblin and D. S. Urch, J. Electron Spectrosc. Relat. 



Phenom. 113, 179 (2001). 

J. P. Rueff, A. Shukla, A. Kaprolat, M. Krisch, M. Loren- 
zen, F. Sette, and R. Verbeni, Phys. Rev. B 63, 132409 
(2001). 

S. Yamaguchi, Y. Okimoto, and Y. Tokura, Phys. Rev. B 
55, R8666 (1997). 

Y. Kobayashi, N. Fujiwara, S. Murata, K. Asai, and H. Ya- 
suoka, Phys. Rev. B 62, 410 (2000). 

R. A. Bari and J. Sivardiere, Phys. Rev. B 5, 4466 (1972). 
S. St0len, F. Gr0nvold, H. Brinks, T. Atake, and H. Mori, 
Phys. Rev. B 55, 14103 (1997). 

For this, it should be confirmed that the IADs and the 
line shape changes are correlated: it will be shown that 
the IAD values of linear combinations of any two normal- 
ized functions are proportional to the IAD value of the 
two functions, and the proportionality factor is the weight 
in the linear combination. Let us consider the functions 
f(x) and g(x), which are normalized to give unit area at 
integration. Consequently, the integral of their difference 
(/ — g) is zero, as integration is a linear operation. The 
corresponding IAD value (let us note it a), however, is dif- 
ferent from zero if / ^ g; thus a = J \f(x) — g(x)\dx. A 
spectrum of an intermediate state can be written as a su- 
perposition of / and g: SP = c\f -\-c2g, where ci +C2 = 1. 
Its difference from the reference g is SP — g = ci(/ — g). 
The integral of the absolute value of the latter is IAD = 
ci J \f(x) — g(x)\dx = cia; this proportionality to the full 
area difference is what makes the IAD value a good mea- 
sure of the extent of the transformation. 



